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ABSTRACT 


This  thesis  presents  the  theoretical  modeling  of  a  terahertz  (THz)  sensor  based  on 
a  flexible  bi-material  cantilever  that  performs  thermal  deformation  under  THz 
illumination.  This  thermal  deformation  was  induced  by  the  large  mismatch  thermal 
expansion  coefficient  of  two  materials.  As  THz  radiation  illuminated  on  the  absorber 
plate,  the  temperature  of  the  structure  increased,  causing  the  cantilever  to  deflect.  The 
magnitude  of  the  deflection  can  be  measured  electronically  or  optically,  which  gives  a 
measurement  of  the  input  THz  signal.  The  bi-material  structure  in  this  thesis  is  composed 
of  a  silicon  dioxide  cantilever  beam,  with  regions  coated  with  aluminum.  The  effect  of 
thickness  of  two  layers  on  the  bending,  response  time,  signal-to-noise  ratio  will  be 
discussed.  Maximum  bending  angle  can  be  found  in  the  results  from  the  example  of 
silicon  dioxide  with  thickness  layer  from  700nm  to  1400nm.  The  imbalanced  stress  due 
to  fabrication  process  was  taken  into  account  to  compare  with  simulation  results  to 
provide  a  better  assumption  on  thickness  selection. 
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I. 


INTRODUCTION 


A.  BACKGROUND 

In  recent  years,  research  has  been  rapidly  developing  in  the  terahertz  (THz) 
frequency  range.  Due  to  its  non-invasive  properties,  terahertz  radiation  has  great  potential 
in  several  areas  of  application  such  as  inspection  of  fragile  artwork,  identification  of 
concealed  weapons,  quality  control  in  manufacturing,  and  medical  imaging. 

Detection  of  THz  radiation  has  been  investigated  and  reported  in  scientific 
literature  [1].  The  sensing  of  THz  radiation  by  the  photonic  detector  techniques  normally 
requires  a  cryogenic  operating  temperature  to  reduce  thermal  noise.  The  imaging  system 
includes  a  vacuum  jacket  to  produce  low  pressure  and  a  cooling  system,  making  it 
relatively  expensive. 

To  overcome  this  problem,  the  potential  use  of  un-cooled  infrared 
microbolometer  camera  technology  for  the  THz  detection  has  been  investigated  [2].  The 
camera’s  original  focusing  optics  were  replaced  with  a  Tsurupica-based  lens  and  the 
milliwatt-scale  quantum  cascade  laser  (QCL)  was  used  as  an  illuminator.  With  this 
change,  the  camera  was  able  to  image  frequencies  down  to  0.7  THz.  One  of  the 
drawbacks  of  microbolometer  technology  is  the  generation  of  heat  (self-heating)  during 
readout,  which  limits  the  integration  time  to  improve  the  signal- to-noise  ratio  [2].  In 
addition,  microbolometer  based  focal  plane  arrays  require  integrated  electronics  for 
readout,  making  them  relatively  difficult  to  fabricate. 

An  alternative  approach  that  eliminates  self-heating  as  well  as  integrated  readout 
electronics  is  detection  based  on  the  deflection  of  bi-material  cantilever  structures  in 
response  to  thermal  radiation  [3,  4].  The  amount  of  deflection  can  be  readout  using  an 
external  light  source.  These  sensors  can  be  easily  fabricated  using  micro¬ 
electromechanical  systems  (MEMS)  fabrication  technologies.  The  bi-material 
micromechanical  sensor  under  consideration  is  based  upon  a  deformed  micro-structure 
cantilever  under  THz  illumination.  The  deformation  is  caused  by  different  rates  of 
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expansion  due  to  the  difference  in  thermal  expansion  coefficients  of  the  two  materials. 
Compared  to  the  photon  detectors,  the  bi-material  un-cooled  thermal  detector  has  lower 
cost  and  is  much  more  portable.  However,  the  relatively  long  thermal  time  constant 
associated  with  high  thermal  isolation  needed  slow  results  in  long  response  times.  It  is 
estimated  that  for  real  time  imaging  (30  frames  per  second),  the  thermal  time  constant 
needs  to  be  less  than  15  ms.  The  main  focus  of  this  thesis  is  to  optimize  the  layer 
thickness  of  thin  films  used  for  fabricating  the  cantilever  structure  to  achieve  high 
sensitivity. 


B.  BI-MATERIAL  CANTILEVER  SENSOR 


The  bi-material  cantilever  deforms  by  absorbing  THz  radiation  by  the  absorption 


plate  (see  Figure  1),  thereby  causing  a  temperature  increase  in  both  materials.  Since  the 
top  material  has  a  higher  thermal  expansion  coefficient  than  the  bottom  material,  the 
cantilever  tends  to  deform  as  depicted  in  Figure  1  [3].  The  magnitude  of  the  deflection  is 
then  measured  electrically  or  optically  and  gives  a  measurement  of  the  input  signal. 


Figure  1.  Deformation  of  bi-material  cantilever  (After:  [3]) 


The  readout  system  can  be  either  electronic  (capacitive  or  resistive),  or  optical. 
Electronic  readout  has  been  the  most  commonly  used  approach,  mainly  in 
microbolometers.  However,  as  mentioned  earlier  it  has  several  disadvantages.  The  main 
problem  is  its  limited  integration  time  due  to  self-heating.  The  electronic  readout  system 
with  a  large  number  of  pixels  requires  complex  circuitry  that  will  make  the 
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microfabrication  process  very  complex  and  expensive.  It  will  also  make  it  very  difficult 
to  scale  the  arrays  to  high  resolution.  Compared  to  the  electrical  readout  system,  the 
optical  readout  is  completely  external  and  does  not  need  any  electrical  leads  attached  to 
the  sensor.  This  allows  the  use  of  low  thermal  conductivity  material,  which  provides 
higher  thermal  isolation  for  achieving  greater  sensitivity.  With  this  advantage,  the  optical 
readout  becomes  more  desirable  for  bi-material  based  focal  plane  arrays  (FPAs). 

Figure  2  shows  the  schematic  diagram  of  an  optical  readout  system  [4].  The 
incident  THz  radiation  from  the  target  source  illuminates  the  backside  of  the  FPAs 
housed  in  a  vacuum  chamber.  The  absorbed  radiation  creates  an  increase  in  temperature 
resulting  in  a  subsequent  deformation  of  the  micro-cantilever  pixels  in  the  array.  On  the 
other  hand,  a  visible  light  source  from  an  LED  illuminates  the  front  side  of  the  FPA.  The 
light  reflected  from  the  FPA  creates  a  thermal  image  and  is  recorded  by  a  CCD  camera. 


THz  soyrce 


Blad€ 


Wrapped 
In  Plastic 


Figure  2.  Schematic  diagram  of  an  optical  readout  system  (From:  [4]) 

As  previously  discussed,  the  main  physical  property  of  the  bi-material  cantilever 
is  a  large  difference  in  thermal  expansion  coefficients  to  obtain  maximum  deflection  and 
achieve  the  operational  goal.  Clearly,  the  absorbing  material  needs  to  have  the  proper 
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absorption  coefficient  to  ensure  the  desired  wavelengths  are  captured.  The  physical 
properties  of  some  of  the  most  commonly  used  materials  available  in  MEMS  technology 
for  bi-material  thermal  detectors  are  shown  in  Table  1  [3,  5,  6]. 


Table  1.  Properties  of  MEMS  compatible  materials  suitable  for  bi-material  sensor 

fabrication 


Thermal  expansion 
coefficient 

a{x\Q-^K-') 

Thermal  conductivity 

g{Wm~'K-') 

Heat 

capacity 

c{Jkg-^K-^ 

Emissivity 

s 

) 

SiNx 

2.1 

19 

691 

0.8 

Si02 

0.4 

1.4 

730 

0.8 

Au 

14.2 

296 

129 

0.02 

Al 

25 

237 

900 

0.02 

As  listed  in  Table  1,  Si02  and  A1  have  the  lowest  and  highest  thermal  expansion 
coefficients,  respectively.  Comparing  physical  properties,  these  materials  provide  the 
largest  difference  in  thermal  expansion  coefficients  and  hence  provide  maximum 
deflection  of  the  cantilevers.  Having  a  low  thermal  conductivity,  Si02  provides  the 
highest  thermal  isolation,  which  makes  it  perfect  for  use  as  the  structural  material  for  the 
sensor. 

In  the  past,  there  have  been  difficulties  in  using  Si02  as  a  MEMS  structural  layer 
due  to  the  lack  of  a  proper  sacrificial  layer  in  the  micro-fabrication  process  [3].  As  an 
alternative,  SiNx  was  used  in  combination  with  Au  or  A1  and  became  a  more  popular 
choice  of  materials  for  bi -material  cantilever  structures  [3]. 

The  Si02  based  thermal  detectors  can  now  be  produced  using  a  new  fabrication 
process  without  the  use  of  a  sacrificial  layer,  as  reported  by  Dong  et  al  [7]  and  as 
discussed  in  a  previous  thesis  work  [3].  Compared  to  SiNx,  overall  deflection 
improvement  using  Si02  was  predicted  to  be  70%  [3].  However,  the  high  stress  of  Si02 
creates  an  extra  curvature  that  affects  the  reflected  LED  light  signal.  This  additional 
curvature  will  remain  a  challenge  but  will  be  discussed  in  more  detail  in  Chapter  IV. 
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II.  SENSOR  DESIGN  AND  MODELING 


The  structure  of  the  micro-cantilever  in  this  thesis  was  based  on  a  previous  design 
that  has  been  widely  used  for  infrared  detectors  [8].  The  two  symmetrical  pairs  of  folded 
legs  hold  the  absorber,  which  also  serves  as  the  reflector  for  the  readout  light  beam,  and 
is  located  at  the  center  of  the  detector  as  shown  in  Figure  3  [6],  Such  configuration 
facilitates  the  magnification  of  individual  deformations.  The  total  deflection  angle  of  the 
reflector  can  be  calculated  as  the  sum  of  the  individual  angles  produced  on  the  folded 
legs,  i.e.  0  =  0^+02- 


Figure  3.  The  geometry  of  micro-cantilever  with  deflection  amplification  (After:  [6]) 


In  order  to  avoid  diffraction  effects,  the  size  of  the  individual  detectors  (pixel)  has 
to  be  larger  than  the  wavelength  of  incoming  radiation  that  is  to  be  detected.  In  this  thesis 
design,  the  measurement  of  one  pixel  is  200//wx200//w  ,  which  is  capable  of  diffracting 
limited  detection  incident  light  from  2  THz  (150//w)  and  upward.  Another  consideration 
of  the  structure  design  was  the  problems  associated  with  partial  device  release  due  to  the 
imperfections  of  the  etching  process.  The  release  involves  removal  of  the  substrate  from 
the  back  side  using  the  Bosch  etching  process  [9].  The  Bosch  process  consists  of  two 
steps:  polymer  deposition  and  reactive  etching  (RIE).  As  shown  below  in  Figure  4  [9], 
during  the  RIE  step,  silicon  was  removed  rapidly  by  the  ion  bombardment  gas.  To 
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prevent  extensive  undercut  during  RIE  (see  Figure  4b),  repeated  polymer  deposition  steps 
are  included  to  protect  the  sidewalls  as  depicted  in  Figure  4. 


Figure  4.  The  steps  of  Bosch  process  where  the  undercut  is  minimized  using  a  thin 

polymer  layer  (From:  [9]) 


This  process  can  result  in  discrepancies  between  the  shape  of  the  mask  of  the 
opening  on  one  end  of  the  wafer  (see  Figure  5a)  and  the  actual  opening  of  the  mask  after 
it  reaches  the  opposite  end,  as  illustrated  in  Figure  5b.  Therefore,  the  design  had  to  be 
adjusted  to  account  for  imperfections  in  the  process. 
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Figure  5.  The  over  and  under  etch  of  different  sections  during  substrate  removal  to 

isolate  pixel 


To  overcome  this  problem,  the  pixel  geometry  had  to  be  redesigned  while  taking 
into  consideration  the  benefits  of  over  etching.  In  this  thesis,  instead  of  having  the  same 
length  of  legs  in  a  square  area  as  in  Figure  6a,  a  shorter  and  outward  pair  of  legs  outside 
the  pixel  was  used  to  preserve  extra  space  for  over  cutting  as  shown  in  Figure  6b.  Figure 
7  shows  the  detailed  dimensions  of  the  detector  designed  by  Dr.  Dragoslav  Grbovic. 


Figure  6.  (a)  Top  view  of  initial  design  and  (b)  top  view  of  new  design  to  reduce  over 

etching  complete  release 
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Figure  7.  Dimensions  of  the  new  pixel  design 


A.  PERFORMANCE  ANALYSIS 


The  goal  of  this  research  is  to  improve  the  performance  of  the  bi-material  pixel. 
The  performance  is  described  by  several  parameters.  The  first  parameter  is  the  per-Kelvin 
deflection  angle  of  cantilever  where  larger  deflection  angle  corresponds  to  higher 
sensitivity.  Therefore,  a  better  performance  can  be  expected  due  to  the  signal’s  being 
stronger  than  the  surrounding  noise.  The  second  parameter  is  the  thermal  time  constant. 
This  describes  how  fast  the  detector  is.  The  heated  material  must  cool  down  to  have  the 
cantilever  return  to  its  original  position  and  be  ready  for  the  next  operational  cycle.  This 
process  will  take  some  amount  of  time,  which  will  limit  the  speed  at  which  the  detector 
can  be  operated.  The  smaller  the  relaxation  time,  the  faster  detector  operation  can  be. 
Therefore,  the  goal  of  this  research  is  to  maximize  the  per-Kelvin  deflection  angle  and 
minimize  the  thermal  time  constant  (response  time). 

To  maximize  the  bending  angle  due  to  the  given  incoming  power,  the  pixel  needs 
to  have  a  large  increase  in  temperature,  AT,  as  given  by  the  following  equation  [3]: 


(2.1) 
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Where  P  is  total  absorbed  power  (watts)  and  G  is  thermal  conduetance  ( W  •  K“' ), 
which  is  the  ability  of  material  to  conduct  heat.  For  the  same  amount  of  absorbed  power, 
the  lower  the  value  of  G  (value  and  based  on  the  material  and  geometry),  the  higher  the 
temperature  change  and  hence  larger  the  deflection  angle. 

Thermal  conductance  also  has  a  great  effect  on  the  heat  conduction.  It  is  inversely 
proportional  to  the  thermal  response  time.  For  the  selected  material  with  a  heat  capacity 
C  {J  the  higher  the  thermal  conductance  G,  the  smaller  (faster)  the  device’s  heat 

dispersion  time.  The  response  time  r  is  determined  by  [3]: 


As  discussed  above,  the  performance  of  the  bi-material  device  is  highly  dependent 
on  the  thermal  conductance  of  the  device.  The  thermal  conductance  depends  on  the 
structure  geometry  design  and  the  materials  used.  The  major  challenge  in  the  design  of 
the  detector  is  the  simultaneous  optimization  of  both  the  bending  angle  and  the  response 
time.  In  order  to  maintain  simplicity  of  processing,  previous  cantilever  based  infrared 
sensor  research  [3]  provided  a  starting  point  for  the  present  work.  The  detailed 
dimensions  of  the  structure  are  listed  in  Table  2. 


Table  2.  Structure  dimensions  used  in  the  analysis 

Pixel  size(//m^) 

Absorber  size 

width  of  legs  ( fxm )  length  of  legs  ( /urn ) 

200x200 

88x200 

8  see  Figure  7 
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B. 


THERMAL  CONDUCTANCE  AND  ITS  LIMIT 


The  equation  for  the  total  thermal  conductance  for  a  bi-material  thermal  detector 
is  given  by  [10]: 


G  = 


+  G  •  7  ^ 

air,conv  air,cond  > 


(2.3) 


Since  the  device  is  typically  operated  under  low  pressure  in  a  vacuum  cell,  the 
convective  (  )  and  conductive  (  )  thermal  conductance  due  to  air  are 

negligible.  The  total  thermal  conductance  can  therefore  be  expressed  as: 

G  =  +  (2.4) 

Gleg  is  the  thermal  conductance  through  the  legs  to  the  heat  sink  and  Grad  is  the 
thermal  conductance  due  to  radiation  emitted  from  the  surface  of  the  pixel.  In  general, 
Grad  is  relatively  small  compared  to  Gieg .  The  value  of  Grad  is  given  by  [1 1]: 

G^^,=4Aar(s„+s„^^)  (2.5) 


where  A  is  the  total  surface  area  of  the  detector,  a  is  the  Stephan-Boltzmann  constant,  T 
is  the  working  temperature  of  the  detector,  and  and  £•5.^  are  the  emissivities  of  A1 

and  Si02,  respectively. 

The  general  equation  for  Gieg  is  given  by  [3]: 


^feg  = 


g-w-t 

I 


(2.6) 


where  g  is  the  thermal  conductivity  of  the  material ,  W  is  the  width  of  the  leg,  t  is  the 
total  thickness  of  the  leg  and  I  is  the  total  length  of  the  legs. 

The  equation  for  the  legs  composed  of  A1  and  SiOi  layers  is  thus  given  by: 


^Al+SiOj 


SiO. 


I 


Al+SiO. 


I 


Al+SiO. 


(2.7) 
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Since  the  is  much  bigger  thang^.^^  ,  the  second  term  can  be  ignored.  The 
equation  becomes: 


Al 


^Al+SiOj 


I 


Al+SiO, 


(2.8) 


The  equation  for  the  legs  with  only  Si02  layer  is  given  by: 


^  _  §Si02  '  ^  '  ^Si02 

^SiO,  ~  ■ 


1. 


SiO, 


(2.9) 


The  for  combining  two  parallel  components  in  2.8  and  2.9  can  be  expressed 
as: 


1  1 

-  +  - 


^leg  ^Al+Si02  ^Si02 


(2.10) 


The  total  for  symmetric  legs  on  both  sides  will  be: 


G,,=2W{ 


S  Al^  AlS  Si02^ . 


SiOj 


^Al+Si02  Ssi02 Gi02  hi02S Al^ Al 


) 


(2.11) 


Again,  since  the  g^,  is  much  bigger  thang^.^^ ,  ignoring  the  term, 

Equation  2.1 1  becomes: 

^Ssi02  ^SiO. 


Gj  ^  2  X 

leg  j 

I  i 


(2.12) 


SiO^ 


Substituting  the  values  of  g^^^^  =lA(Wm  W  =  8//m,  =  282//m 

from  Table  1  and  Table  2  into  Equation  2.7  will  result  in  G,^^  =  3.97  x  10"^t,gg(lEA"') . 

From  Equation  2.5,  the  Grad  can  be  calculated  by  using  ^4  =  1.76x10"* (m^) , 
assuming  the  thermal  radiation  is  only  emitted  from  the  absorber  while  the  surface  area  of 
the  legs  is  too  small  to  be  negligible,  (7  =  5.61  {W~^ K~‘^),  T  =  300J^ ,  and 

is  0.02  and  0.08  respectively  gives  G^^^  =  8.84xlO~^(WK~‘) . 
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Having  a  fixed  value  for  the  width  and  the  length  of  the  legs,  Gieg  ean  be  ehanged 
by  varying  the  leg  thiekness.  The  value  of  Grad  depends  only  on  the  surfaee  area  of  the 
pixel,  whieh  is  nearly  independent  of  the  thiekness  of  the  legs.  Thus,  for  a  fixed  pixel 
area,  the  Grad  remains  eonstant  while  the  Gieg  deereases  with  its  thiekness.  Therefore,  the 
minimum  value  of  thiekness  to  reduee  the  overall  thermal  eonduetanee  ean  be  obtained 
when  Gleg  is  equal  to  Grad  as  illustrated  in  Figure  8. 

However,  due  to  the  issues  of  stress  differential  resulting  from  the  mieofabrieation 
proeess,  700  jum  is  the  lowest  value  of  thiekness  of  Si02  that  has  been  used  in  these 
ealeulations. 
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Figure  8.  Thermal  conductance  as  a  function  of  leg  thickness 
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III.  FINITE  ELEMENT  SIMULATION 


The  response  of  the  pixel  for  absorbed  THz  power  is  simulated  using  COMSOL 
finite  element  simulator.  In  the  model,  the  top  reflector  material  was  temporarily  removed 
in  order  to  eliminate  the  stress  effects  in  the  central  region  of  the  pixel.  In  addition,  the 
exact  material  to  be  used  for  both  absorption  of  THz  radiation  and  reflection  of  the 
readout  is  a  subject  of  a  separate  research  project.  The  value  of  incoming  power  in  this 
simulation  was  estimated  based  on  a  previous  experiment  where  ImW  incident  power 
was  illuminated  on  120x160  pixel  infrared  FPA  with  50 jumx  50 jum  per  pixel  area.  This 

corresponds  to  an  intensity  of  incident  power  of  about  20.8  Wm~^  and  in  COMSOL 
simulation,  100%  absorption  is  assumed. 

It  is  important  to  determine  the  resonance  frequencies  of  the  entire  structure  since 
it  can  affect  the  speed  of  operation  of  the  pixels  under  time  varying  THz  power.  If  the 
lowest  resonance  frequency  is  close  to  the  frequency  of  the  time  varying  THz  signal,  a 
non-linear  response  can  result  from  the  sensor.  On  the  other  hand,  higher  resonance 
frequencies  are  associated  with  stiffer  spring  constant  resulting  in  smaller  deflection  of 
the  cantilever.  Figure  9  shows  the  first  three  eigen  modes  simulated  using  COMSOL  of 
the  pixel  structure  shown  in  Figure  7.  The  first  mode  is  the  most  important  in  the  design 
of  the  pixel. 
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The  frequency  response  of  a  resonance  is  usually  not  be  exactly  symmetric  about 
the  resonant  frequency.  The  amplitude  of  the  oscillations  x  for  a  damped  oscillator  can  be 
expressed  as: 


where  F  is  the  driving  force,  m  is  the  mass,  a>  is  the  driving  frequency,  is  the 
resonant  frequency  and  /  is  the  ratio  of  damping  coefficient  b  to  the  mass  m.  A  typical 
resonant  curve  is  shown  in  Figure  10  (a).  When  the  resonant  frequency  is  much  larger 
than  the  driving  frequency  (i.e.,  co<^cOq),  the  equation  can  be  simplified  as: 

F I  m  F I  m  F 

— —  = - =  —  (3.2) 

k  /  m  k 


which  shows  that  the  response  amplitude  is  proportional  to  the  driving  force,  F  =  kx , 
where  k  is  the  spring  constant.  Under  this  condition,  the  amplitude  of  oscillations  is 
independent  of  the  driving  frequency  and  proportional  to  the  driving  force  providing  a 
linear  response. 


In  the  design,  the  first  mode  for  different  leg  thicknesses  was  simulated  and 
graphed  in  Figure  10(b).  It  was  found  that  the  lowest  resonant  frequency  was  2000  Hz. 
This  result  is  far  from  the  speed  of  operation  of  the  pixels  (30  Hz)  needed  for  real-time 
imaging.  It  is  known  that  a  frequency  of  30  Hz  is  the  minimum  speed  for  a  human  eye  to 
decipher  a  motion  picture.  The  thickness  in  the  simulation  will  present  a  steady  condition 
for  the  amplitude  response. 
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Figure  10.  (a)  Amplitude  as  a  funetion  of  frequency  for  a  damped  oscillator  and  (b) 

resonant  frequency  of  a  pixel  as  a  function  of  leg  thickness 
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A. 


MAXIMUM  BENDING 


As  shown  in  Figure  1 1,  the  total  displacement  of  the  cantilever  in  the  z-direction 
can  be  calculated  by  summing  the  displacement  at  points  1  and  2.  The  sine  of  the 
deflection  angle  (^)  is  trigonometrically  equal  to  the  total  displacement  in  the  vertical 
direction  divided  by  the  length  of  the  pixel.  For  small  0 ,  sin  <9  ~  <9  and  the  deflection 
angle  are  equal  to  the  vertical  displacement. 


Figure  1 1 .  Determination  of  angle  of  deflection 


By  measuring  the  temperature  increase  in  the  absorber  (the  temperature  difference 

dO 

between  heat  sink  and  absorber),  can  be  calculated  for  different  leg  thicknesses. 

Aluminum  has  a  higher  thermal  expansion  coefficient  than  silicon  dioxide  and,  thus, 
increasing  the  thickness  of  the  A1  will  increase  the  bending  angle  as  shown  in  Figure  12 
(a).  However,  it  was  found  that  for  A1  thickness  above  400  nm  (see  Figure  12  (a)),  the 
bending  angle  starts  to  decrease.  This  is  primarily  due  to  reduction  in  temperature  rise  of 
the  pixel  as  the  aluminum  thickness  is  increased  as  a  result  of  higher  thermal  conductance 
of  the  legs.  The  smaller  temperature  increase  incurs  a  smaller  deflection.  The  value  of 
dO 

—  then  can  be  simply  obtained  by  dividing  one  by  the  other,  as  shown  in  Figure  12  (b). 
dT 


16 


dO 

The  —  curve  is  almost  identical  to  the  bending  angle  curve.  As  the  author  fixed  the 
dT 


thickness  of  A1  and  varied  Si02,  the  —  curve  had  the  same  behavior. 

dT 


(a) 


(b) 


Figure  12.  (a)  Deflection  angle  and  change  in  pixel  temperature  as  a  function  of  A1 

dO 

thickness  with  900  nm  Si02,  (b)  dependence  of  —  as  a  function  of  A1  thickness  with 

dT 

900  nm  Si02 
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The  3D  graph  in  Figure  13  shows  the  tendency  of  the  —  curve  by  varying  both 

dT 

aluminum  and  silicon  dioxide  thicknesses.  The  data  in  Figure  13  can  be  used  for 
determining  the  optimum  parameters  for  achieving  the  highest  rate  of  deflection  with 
temperature. 


dO 

Figure  13.  3D  graph  of  the  — by  varying  the  thicknesses  of  both  Si02  and  A1 

dT 


B.  RESPONSE  TIME 

In  order  to  evaluate  the  time  constant  r ,  Equations  2.2  and  2.6  can  be  used  along 
with  the  heat  capacity  of  the  absorber  (200  x  88  pm)  and  thermal  conductance.  By 
assuming  the  mass  of  the  A1  layer  is  too  small  to  be  negligible,  the  heat  capacity  of  the 
structure  can  be  calculated  by  the  specific  heat  capacity  of  Si02  multiplied  by  the  mass  of 
Si02 obtained  from  COMSOL,  C  =  730  {JKg~^K~^)x mass  (Kg) ,  which  means  the 


18 


heat  capacity  is  directly  proportional  to  the  total  thickness  of  the  Si02  layer  since  the 
pixel  area  is  a  constant.  As  mentioned  in  Chapter  II,  the  thermal  conductance  is 
proportional  to  the  thickness  of  Si02  and  Al.  As  a  result,  by  fixing  the  thickness  of  the 
Si02  layer  and  increasing  the  thickness  of  the  Al  layer,  the  time  constant  r  decreased  as 
shown  in  Figure  14.  However,  fixing  the  thickness  of  the  Al  layer  and  increasing  the 
thickness  of  the  Si02  layer,  the  time  constant  r  increased.  This  indicates  that  the  heat 
capacity  dominates  the  thermal  conductance  in  optimizing  the  response  time. 


Figure  14.  3D  graph  of  the  response  time  by  varying  the  Al  and  Si02  thicknesses 

C.  QUALITY  FACTOR 

The  quality  factor,  Q  is  necessary  in  the  calculation  of  the  Noise  Equivalent 
Temperature  Difference  (NETD).  The  details  of  the  NETD  and  its  computation  will  be 
discussed  in  the  next  section.  Q-factor  of  a  resonator  is  generally  defined  as  the  ratio  of 
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the  maximum  energy  stored  to  the  energy  lost  in  a  single  oseillation.  Therefore,  deviees 
with  higher  quality  factors  will  have  larger  stored  energy  when  driven  at  the  resonant 
frequency. 

Q  factor  also  represents  the  sharpness  of  the  resonance  curve  shown  in  Figure  16. 
A  high  quality  factor  is  characterized  with  a  high  response  peak  and  a  narrow  bandwidth, 

which  is  the  frequency  difference  between  the  points  of  amplitude  at  (also  known  as 

v2 

3dB)  of  the  peak,  as  shown  in  Figure  15  [12]. 


Figure  15.  Bandwidth  of  resonant  frequency  (From  [12]) 


Quality  factor  can  be  expressed  as: 


energy  stored  per  cycle  _  resonant  frequency 

energy  dissipated  per  cycle  t  j  •  .  1  r  i 

r'  i'  j-  bandwidth  at  -j=  oj  peak 


(3.3) 


The  damping  constant  (Cdamping)  due  to  the  motion  of  a  cantilever  in  air  can  be 
estimated  using  [13]: 


Champing  =  cantUcver  width  ■  Ti^bpcop 


(3.4) 
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where  co  is  the  vibration  frequency,  p  is  the  fluid  density,  and  p  is  the  fluid  viscosity, 
Z)  is  a  constant  number  dependant  on  the  structure  design.  In  air,  the  density 
is  yC>  =  1.025  k%  •  m“^  and  viscosity  is  //  =  1 .87 1  x  1  •  5  . 

Equation  3.4  was  based  on  a  cantilever  with  rectangular  cross-section  moving  in 
the  air  where  the  parameter  b  is  approximately  equal  to  2.  Considering  the  multi-leg 
design  of  the  structure,  a  higher  damping  and  hence  a  higher  value  of  b  is  expected.  The 
increase  in  b  is  due  to  the  increase  in  the  perimeter  of  the  structure  resulting  in  increased 
air  damping  as  it  vibrates.  The  value  of  b  is  usually  estimated  by  measuring  the 
bandwidth  of  the  resonance  peak  experimentally.  Figure  16  shows  the  experimental 
setup  used  for  measuring  amplitude  as  a  function  of  driving  frequency  using  a  laser 
vibrometer.  A  speaker  was  used  as  a  sound  source,  providing  an  acoustic  plane  wave  in 
the  100  to  7000  Hz  frequency  range.  In  the  measurement,  a  single  pixel  of  a  FPAs 
fabricated  is  employed  with  laser  light  from  the  vibrometer,  focused  using  a  microscope 
objective  lens  assembly. 


Figure  16.  Lab  setup  for  measuring  frequency  response  of  a  pixel.  (A)  Speaker  for 
sound  source  with  directional  cone.  (B)  Laser  vibrometer  with  a  microscope  objective 

lens.  (C)  FPAs 
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As  seen  in  Figure  16,  the  measured  frequency  of  a  pixel  is  compared  with  the 
simulation  with  b  =  3  and  without  residual  stress  (green  line),  showing  a  higher  resonant 
frequency  and  a  symmetric  peak.  However,  a  residual  stress  of  W  Pa  on  A1  and 
3xlO*Pfl  on  Si02  was  found  to  exist  by  the  characterization  of  thin  films  performed 
during  the  microfabrication  process.  Inclusion  of  the  stresses  in  the  simulation  resulted  in 
decrease  of  the  resonant  frequency;  the  curve  became  asymmetric  (blue  line  in  Figure  17), 
which  is  closer  to  the  experimental  data  but  with  a  relatively  wide  bandwidth.  The  red 
line  in  Figure  17  is  the  simulated  data  with  b  =  2  and  stresses  that  show  a  close  agreement 
with  the  measurements. 


Figure  17.  Comparison  of  simulated  frequency  response  with  experimental  data 


By  using  Equation  3.4  with  b  =  2,  the  frequency  response  is  simulated  for  a  set  of 
pixels  with  different  layer  thickness  parameters;  the  results  are  shown  in  Figure  18.  The 
value  of  the  Q  factor  can  be  calculated  by  Equation  3.3  and  ranges  from  20  to  45.  These 
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Amplitude  in  z  direction  (m)  Amplitude  in  z  direction  (m) 


values  are  then  used  to  ealeulate  the  Noise  Equivalent  Temperature  Differenee  (NETD) 
performanee  of  the  deteetor  in  the  FPAs  for  the  next  seetion.  Note  that  the  Q  faetor  in 
operating  eondition  (vaeuum)  will  be  higher. 
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Figure  18.  Frequency  responses  of  pixels  with  Al  thickness  from  160  run  to  200  run 
and  SiOi  thickness  from  800  nm  to  1400  nm 
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D.  NOISE  EQUIVALENT  TEMPERATURE  DIFFERENCE 


Noise  Equivalent  Temperature  Difference  (NETD)  of  a  sensor  is  defined  as  the 
amount  of  temperature  of  a  target  above  the  background  needed  to  produce  a  signal-to- 
noise  ratio  equal  to  1  at  the  sensor.  Both  the  target  and  background  are  assumed  to  be 
blackbody  sources.  An  equation  of  NETD  can  be  expressed  as  [3]. 

7n 

NETD^=—  (3.5) 


where  Zn  is  the  amplitude  of  noise  and  aZ  is  the  amplitude  of  signal.  For  the  micro¬ 
cantilever,  the  signal  will  be  the  deflection  of  the  device  as  a  result  of  power  from  the 
target,  which  can  be  expressed  as  [3]: 

A  rIP 

aZ  =  t,^R(w)(—)^_^T  (3.6) 

where  A  is  the  surface  area  of  a  pixel,  F  is  the  f-number  of  the  lens  used  for  focusing  the 
radiation  beam  onto  the  detector,  and  is  the  transmission  of  the  lens.  In  the  case  of  THz 

detection,  the  lens  is  made  of  polyethylene.  The  average  transmission  of  the  lens  in  the 
THz  region  (1-5  THz)  is  about  0.77,  as  shown  in  Figure  19  [14]. 
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Figure  19.  Transmission  of  2  mm-thick  high  density  polyethylene  in  THz  region  (From 

[14]) 


dP 

R{w)  is  the  responsivity  of  the  sensor,  T  is  the  room  temperature  and  ( — ,  is 

dT  ' 


the  slope  of  the  blaek  body  radiation  curve  ffom/lj  to  .  For  the  desired  frequency  range, 
dP 

1-5  THz,  ( — can  be  obtained  by  approximating  the  power  emitted  by  a  black  body 
dT  '  ^ 

source  in  THz  spectral  range  as  [2]: 


InkT 


vdv 


InkT 


3c 


2  c 


(3.7) 


where  Vc  is  the  upper  limit  of  frequency  range  and  the  lower  limit  is  assumed  to  be  0.  By 
taking  the  derivative  with  respect  to  T  on  both  sides,  the  equation  becomes: 


,dP .  2nk  3 


(3.8) 
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Using  cutoff  frequency  =  5THz  and  calculating ( — )  from  Equation  3.8  gives 

dT 

a  value  of  0.06. 

The  primary  source  of  micro-eantilever  noise  eomes  from  the  temperature 
fluetuations.  The  root  mean  square  (average)  of  this  noise  ean  be  expressed  as  [3]: 

<  Zl,  >‘/2= 

V 


Where  is  the  Boltzmann  constant  of  1.38x10  ^^m^kgs  ,  B  is  the  frequeney 

bandwidth  of  the  sensor — ^which  is  limited  by  the  thermal  time  eonstant  (~  30//z ) — G  is 
the  thermal  conduetance  of  the  sensor  and  r]  is  the  absorbance  (-0.3)  for  a  typieal  THz 
absorbing  material  used  in  sensor  fabrication. 

Additional  noise  is  introduced  from  the  lattiee  vibrations  due  to  the  thermal 
energy  {ksT),  which  excite  some  of  the  resonant  modes.  This  additional  vibration  is  ealled 
thermo-mechanieal  noise  and  is  expressed  as  [3,  15]: 


<Z2  >*'2  = 

^  ^TM  ^ 


\AkJB 


COc. 


Qk 


{col-akf  + 


CO  COa 


(3.10) 


Notiee  that  the  operating  frequeney  co  in  this  case  is  much  smaller  than  the 
resonant  frequeney  ;  therefore  Equation  (3.10)  ean  be  simplified  as: 


^  72  ^1/2 _  l^kgTB 

^  ^TM  ^  ~  , 


Qko}( 


(3.11) 


0 


where  Q  is  the  quality  factor  and  k  is  the  spring  eonstant. 
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By  the  definition  of  NETD  (Equation  3.5),  the  NETD^p  (due  to  temperature 
fluctuations)  and  the  NETD^^  (due  to  thermal  mechanical  noise)  then  can  be  obtained 
by  combining  Equation  (3.9)  with  (3.10)  and  Equation  (3.9)  with  (3.11): 


NETDpp 


AF^T^AKpBG 


(3.12) 


NETDp^ 


4F" 


j4kJB 
V  QkcOo 


(3.13) 


The  total  value  of  NETD  thus  can  be  evaluated  by  taking  a  square  root  of  the  sum 
of  the  squares  of  individual  NETDs. 

dO 

The  graphs  of  NETD  with  —  and  NETD  with  response  time  are  shown  in 

dT 

Figures  20  (a)  and  (b),  respectively.  As  displayed  in  the  graph,  when  the  thickness  of  A1 
is  increased,  the  value  of  NETD  increases  slowly  at  about  0.02  mK  per  50  nm.  When  the 
thickness  of  SiOi  increased,  the  value  of  NETD  increased  rapidly  at  about  2  mK  per  50 
nm. 
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Response  time  (s)  dO/dT  (rad/K) 


Thickness  of  SiO^  (nnn) 
(a) 


(b) 


Figure  20.  (a)  NETD  with  —  (b)  NEED  with  response  time 

dT 
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NETD  (mK)  NETD  (mk) 


IV.  ANALYSIS  OF  RESULTS 


A.  SELECTION  OF  THICKNESS 

It  is  important  to  select  the  most  favorable  thickness  combination  of  the  bi¬ 
material  detector  in  order  to  achieve  the  optimal  performance.  The  challenging  part 
comes  from  the  need  to  simultaneously  optimize  all  the  desirable  parameters  such  as 
sensitivity,  response  time  and  noise.  As  shown  in  the  equations  above,  these  complex 
physical  properties  have  adverse  effects  on  each  other.  Therefore,  optimizing  one  often 
worsens  another. 

One  way  to  simplify  this  problem  is  to  create  a  simulation  of  numerous  values  of 
different  thickness  settings  and  then  graph  the  results.  The  desired  properties  of 
unsimulated  data  can  be  predicted  from  the  resultant  curve.  From  previous  chapters,  the 
3D  graph  in  Figure  14  shows  that  the  response  time  will  decrease  by  either  reducing  the 

dO 

thickness  of  Si02  or  increasing  the  thickness  of  Al.  The  maximum  value  of  —  can  be 

dT 

found  in  each  curve  from  the  3D  graph  on  Figure  13,  which  was  simulated  by  100  nm 
intervals  of  the  Si02  layer  thickness.  The  peak  point  with  the  highest  deflection 
amplitude  and  the  smallest  total  thickness  with  the  smallest  response  time  from  Figure  21 
will  provide  the  best  solution.  In  conjunction  with  these  selection  criteria,  the  least  total 
thickness  that  can  be  considered  must  be  700  nm  due  to  restrictions  on  current  fabrication 
techniques.  As  discussed  previously,  the  value  of  NETD  of  all  the  combinations  with 
Si02  thicknesses  less  than  1400nm  will  never  be  higher  than  62  mk,  which  is  relatively 
low  for  uncooled  thermal  sensors  in  the  THz  spectrum. 
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Figure  21 .  Response  time  and  peak  defleetion  rate  as  a  funetion  of  Si02  thiekness 


Another  method  for  determining  the  eombination  of  thieknesses  is  to  fix  one  or 
two  variables  and  vary  the  remaining  parameters.  For  example,  in  Figure  22,  if  the 
desired  NETD  must  be  lower  than  that  of  the  thiekness  eombination  of  4  (Si02  800  nm, 
A1  200  nm),  whieh  is  46  mK,  and  the  response  time  has  to  be  lower  than  thiekness 
eombination  1  (Si02  800nm,  A1  80nm),  whieh  has  a  value  of  0.75  seeond,  this  will  limit 
the  ehoiee  of  bending  angles  of  eombinations  2  and  3.  The  eombination  of  3  will  thus  be 
seleeted  beeause  it  provides  the  maximum  angle. 
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Figure  22.  Thickness  sets  with  the  optimized  desirable  parameters 


B.  THE  EFFECT  OF  STRESS  FORCE 

Another  important  physical  effect  on  the  performance  of  bi-material  cantilevers  is 
the  initial  deflection  due  to  the  residual  stress  in  the  two  materials.  This  stress  builds  up 
during  the  chemical  deposition  stage  of  the  micro-fabrication  process  [16,  17].  The 
residual  stress  is  a  result  of  a  difference  in  the  shrinking  rate  between  two  materials  in  the 
x-direction  and  y-direction  when  cooled  down.  These  effects  form  a  curvature  as  shown 
in  Figure  23. 
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(a) 


(b) 


Figure  23.  Top  and  side  views  of  (a)  structure  without  residual  stress  and  (b)  structure 

with  residual  stress 


The  initial  curvature  highly  affects  the  optical  readout  caused  by  thermal 
deflection  due  to  THz  absorption.  These  factors  are  important  for  consideration  in 
performance  analysis  in  order  to  properly  choose  the  realistic  thicknesses  for  simulation. 
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V.  CONCLUSION 


This  thesis  has  successfully  demonstrated  the  use  of  COMSOL  multi  physics 
finite  element  simulation  with  heat  transfer  and  solid  stress  models  for  the  design  of  a  bi¬ 
material  cantilever  sensor  in  the  THz  region.  The  results  for  changing  thicknesses  and  the 
physical  effect  on  the  sensitivity,  response  time  and  noise  were  reported.  From  these 
results,  the  best  combination  will  be  compared  experimentally  to  the  proposed  device. 
The  characterization  methods  developed  in  this  study  laid  a  foundational  process  for  a 
more  advanced  structure  design  and  more  complex  simulations.  In  the  future,  coupling 
these  results  with  the  COMSOL  RF  model  could  enable  a  simulated  heat  generation 
source  from  THz  radiation  to  more  closely  emulate  the  actual  device. 
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